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ABSTRACT 


This report describes results of studies conducted by the Research 
Triangle Institute under Contract NAS1-11298. An analysis of Very Low 
Frequency propagation in the atmosphere in the 10-14 kHz range leads to a 
discussion of some of the more significant causes of phase perturbation. 

The method of generating sky-wave corrections developed by the Naval 
Electronics Laboratory Center to predict the Omega phase is discussed. 
Composite Omega is considered as a means of lane identification and of 
reducing Omega navigation error. A simple technique for generating 
trapezoidal model (T-raodel) phase prediction is presented and compared 
with the Navy predictions and actual phase measurements made by J. A. Pierce 
at Harvard University. The T-model prediction analysis illustrates the 
ability to account for the major phase shift created by the diurnal effects 
on the lower ionosphere. The T-model predictions exhibit promise as a means 
of generating phase corrections due to the inherent simplicity of their 
generation. Finally, an analysis of the Navy sky-wave correction table is 
used to provide information about spatial and temporal correlation of phase 
correction relative to the differential mode of operation. 
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1.0 INTRODUCTION 


The need for simple and precise navigation systems which provide large 
area coverage is rapidly expanding as the number of aircraft and ships 
increases. Some presently-used aircraft navigation aids (e.g., DME) are 
active devices which have limited range and definite saturation limits in 
today's growing population of aircraft. Passive techniques which cannot be 
saturated and are potentially low in cost will have an increasingly important 
role in the future of navigation. The Omega very low frequency (VLF) naviga- 
tion system developed by the Navy as a means of long-range navigation for ships 
at sea offers many advantages as an aircraft navigational aid. World-wide 
coverage can be attained with only eight transmitter stations, and consistent 
accuracies of less than 2 nautical miles (n. mi.) rms have been demonstrated in 
the normal operational mode. In a differential mode, consistent daylight 
accuracies on the order of .5 n. mi. rms are reported for ranges up to 300 n. mi. 

A need exists to determine the limiting factors in VLF navigation 
accuracies and to explore new techniques to use all of the information in the 
Omega format of transmissions to improve accuracies. To fulfill this need, 
more knowledge must be obtained about the propagation characteristics of the 
VLF signals, the factors which contribute to phase perturbations, and the 
anomalies which affect phase stability. Additionally, simple techniques are 
needed to apply the necessary corrections to raw phase and phase difference 
readings so that the system may be available in an inexpensive, simple-to- 
operate version for use by general aviation aircraft in situations where 
navigation aids such as VHF Omnirange (VOR) and Distance Measuring Equipment 
(DME) are not available and Visual Flight Rules (VFR) navigation may not be 
practical. 

This report describes work which has been directed toward an investigation 
of a very simple technique for providing Omega phase measurement corrections. 
Additionally, some effort has been applied toward the concept of reducing 
navigation error by using composite frequency phase measurements. Finally, 
an initial look at the technique termed Differential Omega is considered, 
not only to provide a new perspective on analyzing the medium, but as a means 
of significantly reducing the Omega navigation error. 

Chapter 2 discusses the Omega navigation system and provides some back- 
ground material. A bibliography has been compiled in Appendix A which 
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includes a large number of publications pertaining to VLF propagation and the 
Omega navigation system. Chapter 3 provides a summarized analysis of the 
waveguide model, a description of the known factors contributing to phase 
perturbations, a summary of the method of prediction of sky-wave corrections 
(SWC) used and published by the Navy, and an analysis of some aspects of 
composite frequency phase propagation characteristics. 

Chapter 4 describes the trapezoidal model for generating sky-wave 
corrections, and Chapter 5 provides an analysis of this model, comparing it 
to the Navy predictions and to some actual phase measurement data obtained 
from Professor J. A. Pierce of Harvard University* Finally, Chapter 6 discusses 
the differential Omega concept and provides some insight, through analysis of 
sky-wave correction tables, into factors affecting differential Omega accuracy. 


^Permission to use the Cambridge data was granted by Professor Pierce. 



2.0 THE OMEGA NAVIGATION SYSTEM 


2 . 1 General 

The Omega navigation system is a world-wide hyperbolic system which will 
ultimately employ eight transmitter stations providing a global navigation 
capability on the earth's surface and in the air. Presently, there are four 
transmitters in operation at frequencies of 10.2, 11^, and 13.6 kHz. The 
present stations are located at Aldra, Norway; Trinidad, West Indies; Oahu, 

Hawaii; and in North Dakota. A station at Shashima, Japan, is currently under 
construction and plans have been made to provide transmitters at Reunion Island 
in the Indian Ocean, at Tre-lu, Argentina, and in Australia. Eventually, all 
stations will have output power of 10 kW and will use atomic frequency standards. 

The U. S. Navy has been primarily responsible for the development of the 
Omega system and through the administration of the Omega Projects Office has 
designed the specifications for equipment and the format for transmissions. 
Recently, the U. S. Coast Guard has been designated to administer the Omega 
system and the Defense Mapping Agency is responsible for publishing the sky-wave 
correction tables previously published by the Navy. 

2.2 How Omega Works 

Each transmitter station transmits each frequency for approximately one 
second (continuous wave) every ten seconds according to a prescribed format. 

Table 1 is a chart which shows the eight station formats. The entire format 
is repeated every ten seconds. 

The present Omega configuration, consisting of four transmitting stations 
in the northern hemisphere, provides coverage of one third of the world with a 
theoretical probable positioning error of + 1 n . mi. daytime and + 2 n. mi. night- 
time. Signals from transmitter stations can be received out to maximum usable 
ranges of 4000 n. mi. to 8000 n. mi. from the transmitters, depending on the 
bearing of the receiver from the transmitter. The shorter range can be 
expected west of the transmitter, the longer range east of the transmitter. 

The system is based on the principle that the velocity of propagation of 
electromagnetic radiation is constant. Thus, the difference in distance from 
a receiver to each of two synchronized transmitting stations is indicated by 
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Omega Stations 


Station Name 

Designator 

Latitude 

Longitude 

Norway 

A 

66°25’15.00"N 

13°09 1 10.00"E 

Trinidad 

B 

10°42 , 06.20 M N 

61°38 , 20.30 ,, W 

Hawaii 

C 

21°24 t 20.67"N 

157 0 49’47.75"W 

New York 

D 

43°26’40.92"N 

75°05 , 09.80"W 


Table 1. Omega signal format. 
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the difference in the phase of the signals received. This phase difference 
may be a fraction of a wavelength or many wavelengths, depending upon the 
relative distance between the receiver and each transmitter. The locus of 
points at a constant difference in distance between the receiver and trans- 
mitter stations is a hyperbola. At any given constant difference in distance, . 
a hyperbola or line of position (LOP) is described. All LOP's intersect 
the line joining the transmitter stations (baseline) at right angles. Thus, 
a receiver measuring the difference in distance between two transmitters can 
be placed somewhere on a particular LOP between the transmitters. 

Using a second pair of transmitters, measurement of phase difference 
will place the receiver on another LOP. The intersection of the two LOP's 
then defines a point on the earth's surface which determines the location of 
the receiver. The navigator will choose Omega LOP's for good accuracy and 
for large crossing angles. 

There are an infinite number of hyperbolic lines of position between 
any two transmitting stations. With the transmitter sites fixed at precise 
geographic locations and all transmitting signals accurately phase synchronized, 
the hyperbolic LOP's formed between each pair of stations are 
fixed and precisely located. Figure 1 provides a sample of published 
navigational charts which show a representative quantity of these LOP's. 

These charts are also provided in tabular form (see ref. 1) so that each LOP 
is defined in terras of a set of latitudes and longitudes defining points 
through which each LOP passes. The navigator can then determine from a pair 
of LOP measurements the latitude and longitude of his receiver position. 

Since phase difference measurements are cyclic as the receiver position 
changes, there is in practical application ambiguity in determining the 
position of the receiver. Lanes between LOP's corresponding to the same 
measured phase difference represent the largest regions of unambiguous 
position determination. These lanes have a ground distance width along the 
base line corresponding to one-half wavelength at the particular frequency 
used. (See Appendix B.) For example, the LOP lane width is approximately 
8 n. mi. at 10.2 kHz. Thus, in using Omega, the navigator must know his 
receiver position within 8 n. mi. to obtain an unambiguous position fix from the 
LOP measurements. At the higher frequencies the lane widths are reduced 
somewhat, but by using two frequencies, such as 10.2 and 13.6 kHz, the 
composite frequency 3.4 kHz lane width of unambiguous position determination 
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is extended to 24 n. mi. Consequently, by using the 1.133 kHz difference 
between 10.2 and 11^- kHz, it is possible to extend the lane width to 72 n. rai. 

2.3 Using Omega to Determine Position 

As an example of how Omega phase measurements can be used to determine 
position, consider the following illustration. A.Tracor Omega receiver 59 9R, 
located in Building 1299 at Langley Research Center (LRC) , Hampton, Va., was 
used to make measurements. The exercise illustrates how the published sky- 
wave correction (SWC) tables are used to provide for a position fix. 

On May 11, 1972, at 1300 hr. Omega phase measurements with respect to 
the 599R internally generated reference signal were made at 10.2 kHz. 

Table 2 summarizes the phase measurement readings. 


Table 2. Selected station relative readings, 
(centicycles @ 10.2 kHz) 


Time/Date 

Stations 

A-R 

B-R 

C-R 

D-R 

F-R 

1300/5-11-72 

84 

20 

24.5 

27.5 

87.5 


Six phase differences calculated from the data in Table 2 are given in 
Table 3. 


Table 3. Station pair readings. 


(centicycles @ 10.2 kHz) 


Time/Date 

Stations 

A~B 

A-C 

A-D 

B-C 

B-D 

C-D 

1300/5-11-72 

64 

59.5 

56.5 


gg| 

M ■ 


Using the SWC tables (see ref. 2), the corrections for 1800Z on 11 May 1972 
for the stations are repeated in Table 4. The location of the receiver is 
estimated at 37°5'N, 76°23'Wso that the SWC for 36°N, 76°W are used. 
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Table 4. SWC for 36°N, 76°W. 


Time/Date 

Stations 

A 

B 

C 

D 

1800Z/5-11-72 

-4 

5 

-1 

-7 


Station Pairs 

A-B 

A-C 

A-D 

B-C 

B-D 

C-D 

-9 

-3 

3 

6 

12 

6 


The corrected readings for the station pairs are tabulated in Table 5. 


Table 5. Corrected readings. 


Time /Date 

Station Pairs 

wsssm 

■Kssxai 

K53EI 

■€9301 



180QZ/5-11-72 

LOP 

55 

1001.55 

56.5 

845.56 

59.5 

1088.59 

1.5 

744.01 

4.5 

987.04 

3 

1096.03 


Using the Omega chart (see ref. 3), it appears that the receiver at LRC was 
located in lanes AB 1001, AC 845, AD 1088, BC 744, BD 987 and CD 1096. Using 
the results presented in Table 5, the corresponding LOP’s of the receiver are 
given in Table 5. 

Using the Omega lattice tables, the Omega lines can be plotted on a map. 
For pair B-C, the corrected reading is 744.01. From the B-C lattice tables 
(see ref. 1, p. 48), 

Lat . Tabulated Long. A_ 

37°N 76.28.4W 13.3' 

38°N 75.32.6W 13.8' 

the lane is BC 744. Thus, the corrected reading less the lane count is 

.01 cycle. To find the interpolated longitude, multiply A and the difference 


*Corrected. 
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Of 0.01. Then 


(13.3 ' ) (.01) = . 13’ 
(13. 8' ) (.01) = .14' 


and 


so that the interpolated longitude becomes 


76°28.5’ (? 37°N , and 
75°32 . 7 ' @ 38°N . 


Table 6 summarizes these results and the interpolated LOP points for pairs BC, 
BD, AB, and CD. 


Table 6. Interpolated LOP’s for Omega fix. 


Station 
Pair LOP 

BC 744.01 

BD 987.04 

AB 1001.55 

CD 1096.03 

Long., Lat. 
Points of 
Intersec- 
tion 

76°28.5’W37°N 

76 o W37°10’N 

76°W37°6.7’N 

76°W36°41.1'N 

75°32.7’W38°N 

77°W37°00.8 t N 

77°W37°5 ,8’N 

77°W37°3.5'N 


Figure 2 is an overlay plot to determine the Omega fix of the receiver. The 
LOP’s BC, BD, and AB intercept at a point 37°6'N 76°23’W. Line of position 
CD 1096.03 is several miles off. This discrepancy cannot be explained except 
that the CD pair is not recommended for use in this area. 
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Fig. 2. Plot of the Omega fix. 
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3.0 ANALYSIS OF VLF PROPAGATION 


3 . 1 Waveguide Model 

Determining the electromagnetic field radiated between the earth and 
the ionosphere by an Omega transmitter is basically a boundary-value problem. 
The field is expanded in wave equation solution functions which are nearly 
compatible geometrically with the earth-ionosphere waveguide and then boundary 
conditions are applied to provide equations which may be solved for the eigen- 
values of the solution functions. Unfortunately, the problem is an extremely 
difficult one because the boundary conditions and shape of the guide vary in 
such a way as to make a closed form solution impossible. All available solu- 
tions, therefore, have simplifying assumptions as cornerstones in their 
developments. The idea in solving for the fields radiated by the antenna is 
to determine which assumptions can be made without invalidating the results. 

Although the problem of an antenna radiating into the earth-ionosphere 
waveguide is one of the oldest electromagnetics problems, not until 1959 did 
Budden (ref. 4) offer a method to obtain a full wave solution to a model 
which takes into account a curbed earth, a ground of finite but constant 
conductivity, an ionosphere whose properties vary with height above the earth, 
and an arbitrary but constant orientation of the earth’s magnetic field. 

Budden ’s model was actually rectangular with the earth's curvature being 
taken into account by a varying permittivity. A computer program to calculate 
the fields using this model was developed by Pappert et al. at NELG 
(ref. 5), reported in 1967. Prior to this, in 1964, Wait and Spies (ref. 6) 
offered a solution using a spherically-shaped waveguide model and a layered 
ionosphere which also took into account all the critical parameters. Their 
method of including the effect of the anisotropy of the ionosphere was less 
rigorous, however, and although the first order mode calculations are very 
close between the two methods, there is a significant difference in second 
order calculations. The work of Wait was directed toward obtaining a closed 
form mathematical solution for ease of calculation, whereas the Pappert 
formulation is more dependent on the computer. 

Even though the calculations of Pappert take into account the inhomo- 
geneity and anisotropy of the ionosphere, the waveguide model is nevertheless 
a homogeneous one. This means that the properties do not change along the 
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direction of propagation in the guide which is not true of the ionosphere; 
thus it would be inconsistent to attempt a rigorous full wave solution of the 
entire ear th-iono sphere waveguide. Therefore, assumptions are employed in 
order to use the homogeneous guide results in a practical problem. 
Consequently, the chief use of the calculations and theory based on a homo- 
geneous waveguide has been to assist in the design of more practical 
propagation models such as the one due to E. R. Swanson of NELC, which is 
used to generate the tables of correction values through the U. S. Naval 
Oceanographic office. 


3.2 Factors Which Affect Omega Phase 


To "affect" phase means to alter it, and in order to measure the 
alteration a standard is first needed from which a deviation can be 
noticed. For reasonable interpretation it is believed that the unaltered 
standard should be the first order mode of the Omega wave propagated in 
a spherical earth-ionosphere waveguide with no magnetic field, a nominal 
ionosphere varying exponentially with height only, and a homogeneous ground. 
Under these conditions the phase waveguide problem has been solved, and 
changes in conditions will alter the phase from the predicted value. 
Practical considerations, of course, dictate the types of changes that 
can occur and the ones that will affect the phase appreciably. 

The propagation path from Omega transmitter to receiver is the region 
wherein changes in the earth-ionosphere waveguide will affect the received 


signal. This path is principally the first fresnel zone of the receiver 
looking back at the transmitter. The central width may be approximated by 
(see ref. 7) 


cos 


w 

2a 


COS % + 2l> 


cos 


2a 


where a = radius of earth, 

d = distance between transmitter and receiver, 
and X = Omega wavelength. 

The causes of phase changes can be placed into two categories: one, 
those caused by higher order mode interference; and two, those caused by 
changes in the conditions of propagation of the first mode. The higher order 
mode interference occurs near the transmitter and near boundaries in the 

propagation path. The higher order modes excited at the transmitter attenuate 
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quickly so that they are of little consequence at distances over 650 n. mi. 
from the transmitter (400 n. mi. for all land paths). At boundaries in the 
propagation path such as a sea-land interface or day-night line, first order 
mode waves may couple into higher order modes. Fortunately, mode coupling 
at the lower end of the VLF spectrum, appropriate to Omega, is slight and the 
higher order modes are attenuated quickly and are only of local consequence. 

The first order mode attenuation and phase velocity become the 
factors of interest to a navigating user of the Omega system. Atten- 
uation rates of VLF are very low and are therefore only a factor during 
anomalous propagation conditions. Sunspot activity and solar flares can 
cause absorption of VLF waves at the higher latitudes. An enhancement of 
solar cosmic rays (SCR event) can cause both auroral absorption (AA) and 
polar cap absorption (PCA) . Auroral absorption occurs between 62.5° 
and 65° geomagnetic latitude and can last up to two weeks. Propagation 
from the Norway Omega station at 66°N latitude would be the only station 
affected by this anomaly for users south of the auroral absorption zone. 

Polar cap absorption is also caused by SCR's and affects principally LF , 
and sometimes VLF, waves above 60° geomagnetic latitude. For propagation 
paths below 60®, attenuation should not be a factor, although reception of 
the signal may be altered by variables such as precipitation static and 
interference from noise sources. 

The phase velocity is the limiting factor for Omega accuracy. It is 
affected by variations in ionospheric height which in turn vary seasonally, 
daily, and with latitude, as well as with propagation anomalies caused by 
solar flares. Phase velocity is also affected by ionospheric roughness 
and by the complex reflection occurring at the ionosphere that is a 
function of the relation between lower ionospheric electron density profile 
and the geomagnetic dip angle. One more factor which fortunately seems to 
be independent of the others is ground conductivity. Despite all of these 
factors the phase of the Omega wave is extremely stable, within 6 centicycles 
(cec), and predictable to within a phase difference shift of 10 to 20 cec. 

Ten cec are normally equated to position error of 1 n. mi. when navigating 
using phase differences (see Appendix B) . Changes in phase velocity due to 
solar activity cause the predictability to be much worse, however. Ionospheric 
lowering due to sudden ionospheric disturbances (SID) caused by X-radiation 
from a solar flare can cause errors of up to 30 or 40 cec (3 to 4 n. mi.) for 
all daylight paths. These disturbances form in a matter of minutes after 
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the flare begins and last 2 to 3 hours. Extreme flares ("M" type flares) 
have been recorded that caused errors for 12 hours, the maximum being 
20 cec on a Trinidad to Forrestport, N.Y., path. Although rare, they are of 
significance because of their duration (see ref. 8). PCA events have much 
the same effect on the phase and are of longer duration, 5 to 15 days, but 
they are only of consequence for propagation paths that extend into the 
higher latitudes. For most navigation situations, other lower latitude 
Omega stations could be selected. 

3.3 Naval Electronics Laboratory Center (NELC) 

Phase Prediction Program 

A computer program to predict the phase of the 10.2 kHz Omega signals 
in selected regions of the earth was generated at NELC in order for Omega 
navigators to correct their phase readings for propagation conditions. 

E. R. Swanson described the program as using a single mode parametric model 
with force-fit statistical smoothing (see ref. 9). Only the principal 
waveguide mode is considered, and parameters which describe the propagation 
are directly related to easily defined path characteristics such as orienta- 
tion, latitude, ground conductivity, and time of day and year. 

Electromagnetic waveguide solutions would relate the propagation 
parameters to physical characteristics of the guide such as ionospheric 
height, electron density profile, magnetic field, and ground conductivity. 

The parametric approach was taken because of the uncertainty in these 
characteristics, although the waveguide solutions were used to provide 
guidance on the separability of the variables and to suggest functional 
forms for the variations. The parametric approach amounts to postulating 
functional forms with unset dummy parameters for the variation of the 
phase with the path characteristics. These dummy parameters are then 
chosen by attempting to isolate each one's effect on the phase and 
picking the one which gives the best fit to measured phase data. Regression 
analysis techniques are used to determine the generating function parameters. 

Parameters based on data collected over 11 years from 300 sites, along 
with continuous modification of the phase prediction program, have contributed 
to make the computed position of the navigator using the Omega sky-wave 
correction tables to be within 1 mile of his actual position (ref. 10). 
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3.4 Composite Omega 


Composite Omega refers to the use of multiple frequencies from the trans- 
mitter stations in making phase difference measurements. The ability to take 
advantage of the available phase information at more than one frequency can 
serve to enhance navigation accuracies in that the phase perturbations at VLF 
are to some extent frequency dependent (refs. 11-13). Furthermore, as has been 
mentioned in Chapter 2, the use of difference frequency phase (between pairs 
of transmitter Omega frequencies) can serve to increase the physical width 
of the region of unambiguous phase difference measurements. With lane width 
increases, the utility of Omega navigation increases, since the navigator 
does not need to know his approximate position as accurately as in single-frequency 
Omega. As mentioned in Chapter 2, the lane width of the 3.4 kHz* beat frequency 
is 24 n. mi. as opposed to the 8 n. mi. lane width of the 10.2 kHz signal (see 
Appendix B) . 

The extended lane width accomplished through the use of composite Omega 
makes it possible to use a difference frequency to identify in which 10.2 kHz 
lane the receiver is located. For example, the 3.4 kHz signal phase could, in 
principle, identify one of three possible 10.2 kHz lanes, thus extending the 
interval between ambiguities to 24 n. mi. As noted in ref. 13, the primary 
problem in lane identification is the error in prediction of phase. 

Since transmission time T^ at a given frequency f is proportional to 
phase over a given path, the relationship 


T 3 .4 4 T 13.6 “ 3 T 10. 2 


( 1 ) 


relates the phase of the 3.4 kHz difference frequency to the phase of the 
13.6 kHz and phase of the 10.2 kHz Omega signals (see Appendix C) . In (1) , 
all times or phases are in units of centicycles (cec) of 10.2 kHz, 

If, as an example, an unexpected variation of 9 cec at 10.2 kHz 
occurs and a variation of 5 cec in the same direction and in the 
same units occurs at 13.6 kHz, then from (1), 


AT 3 4 = -7 cec of 10.2 kHz (2) 

occurs at 3.4 kHz. The difference between the 3.4 kHz and 10.2 kHz signals 


*3.4 kHz = 13.6 kHz - 10.2 kHz. 
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3. Omega 10.2 kHz, 11.3 kHz, and 13.6 kHz signals initially in phase and resultant sum and 
quadrature sum signals. 


has changed by 16 cec. Thus, as pointed out in refs. 11 and 13, part of the 
lane identification problem occurs because propagationally-induced variations 
in 10.2 kHz and 3.4 kHz usually appear in opposite senses, and the difference 
is greater than either change alone. For this reason it is particularly 
critical that predictions be very good in order to obtain reliable lane 
identification with composite Omega. 

Langley Research Center has indicated some interest in a three-frequency 
composite Omega system where the sum signal of all three frequencies is used 
and permits resolution of ambiguities over approximately 72 n. mi. on the 
baseline between stations. 

To illustrate the sensitivity of such a technique to phase perturbations 
of individual Omega frequencies, consider a detector which might track the 
sum signal using a phase-locked loop. The zero point of the sum signal where 
the individual component signals are also zero and positive going represents 
the point where the phase detector is locked on, as illustrated in Fig. 3. 
This point reoccurs after 9 cycles of the 10.2 kHz signal, 10 cycles of the 
11.3 kHz signal, and 12 cycles of the 13.6 kHz signal. The quadrature sum 
signal is also shown on the figure. The horizontal scale is in nautical 
miles along a baseline between station pairs, and permits resolution of 
ambiguities over approximately 72 n. mi. on this baseline. 

Figures 4 and 5 show the effect of a phase shift on one or more of the 
individual phase detector outputs. In Fig. 4 the 11.3 kHz signal has been 
shifted -43 deg. and, as may be seen, an error of approximately .2 n. mi. 
appears on the zero point of the sum signal. In Fig. 5, the 13.6 kHz signal 
has been shifted -90 deg. and the sura signal zero point has changed by almost 
one complete cycle of the 10.2 kHz phase, or 7.4 n. mi. 

While the three-frequency composite Omega system is tolerant of small 
uncorrected phase variations, larger uncorrected phase variations of the 
individual signals will cause an automatic phase tracker to, in effect, 
skip cycles and provide erroneous position data. Thus, as has been pointed 
out in two-frequency composite Omega, care must be taken in correcting each 
of the individual frequencies in three-frequency composite Omega. 

As indicated in refs. 11 to 13, composite Omega in which the phase data 
from two separate carriers are linearly combined in a specific manner can 
exhibit higher phase stability and predictability than that obtainable with 
either carrier alone. Also, the composite signal offers potentially greater 
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accuracy than the simple 3.4 kHz difference frequency (ref. 11). In ref. 12, 
conclusions made on the basis of data analysis indicate that use of this 
composite signal can reduce diurnal variation, propagation anomalies, and may 
make it possible to obtain sufficient accuracies without the use of sky-wave 
corrections. This would have a profound effect on the use of Omega in aircraft. 

In ref. 11, the idea of using the composite signal to detect the 
presence and magnitude of propagation anomalies is stated. In a mode when a 
difference frequency is used to make lane identification, the known magnitude 
of an unpredictable phase shift could be used to correct both the carrier and 
difference frequency to nearly the values they would have in the absence of 
an anomaly. Thus, lane count could be maintained while taking advantage of 
the increased precision of the composite signal. 
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4.0 TRAPEZOIDAL CORRECTION MODEL FOR OMEGA (T-MODEL) 


4.1 Introduction 

The trapezoidal model (T-model) (ref. 14) for determining predicted transmitter 
phase measurement at a given receiver site is one of the simplest means of 
providing such information. When compared to the complexity of calculating 
the Navy sky-wave correction tables (see ref. 9) or to the polynomial 
generating function (see ref. 15), the T-model offers a great advantage. 
Fundamentally, the T-model accounts for the diurnal variations by using, an 
all-daylight condition phase prediction, an all-night time condition phase 
prediction and assumes a linear change in phase during transition times. The 
critical parameters in completely describing the T-model for a given 
transmitter-receiver path are the sunset and sunrise times at the receiver 
and transmitter sites, the range between the transmitter and receiver, and 
the relative phase velocity at night and during the day. The relative phase 
velocities used are generally long-term averages determined through experimen- 
tation. 

This chapter describes a simple method for generating the T-model phase 
prediction and sky-wave corrections. The method presented is based in part 
on previously documented schemes but has been simplified and generalized. 

Needed inputs are the latitudes and longitudes of the transmitter and 
receiver, the sun ephemeris, an "average" earth radius, and the estimates of 
daytime and nighttime relative phase velocity of the Omega signal at the 
frequency of interest. Corrections are determined by comparing phase predic- 
tions with the nominal or chart-value* of phase. Finally, a discussion of 
error in T-model predictions as a function of errors in estimating transmitter- 
receiver path lengths is provided. 


4.2 Calculation of Sunset and Sunrise Times 


Given the latitude and longitude of any point on the earth, the sunrise 
and sunset times referenced to Greenwich mean time (GMT) can be calculated. 
To accomplish this, the given latitude (map latitude), which is a geodetic 
latitude, must be converted to a geocentric latitude which, geometrically 


*Chart-value refers to phase computed using (c/v , ) = .9974 as in 

the Navy SWC charts. P 36 
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speaking, describes the given latitude on the earth's surface as a projected 
latitude on a spherical approximation to the earth's surface. 

The geodetic system describes the latitude of a point on the elliptical 
earth's surface in terms of an angle measured between the major earth axis 
and an intersecting line drawn perpendicular to the earth's surface tangent 
line at the point of interest. In Figure 6, <j> represents the geodetic latitude 
of point (x,y) where a cross-section view of the earth is illustrated. From 
the equation of an ellipse, 



To find the slope of the tangent to the ellipse at point (x^y^, take the 
derivative such that 


or 


2xdx 

2 

a G 


+ 


2 2*L 



= 0 


d£ 

dx 
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x 

y 
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Earth ' s 
Surface 


Elliptical Earth 


Fig. 6. Cross-section of the earth's surface showing geodetic and 
geocentric latitudes . 
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But 


i2L 

dx 


-cot <J> and — = cot 0 * 


so that 


0 

tan 4 > = — — tan <j> ' 


or 


* _ 


= tan 


a[V 

L a n" 


tan $ 


( 1 ) 


Equation (1) describes the geocentric latitude of point (x^y^) , 4 ’ , in 

terms of the geodetic latitude <J> and the square of the ratio of minor to 

major axis of the earth. Using the values of and given in ref. 16 which 

are based on the Clark Spheroid of 1866 (used as the basis for U.S. geo- 

2 2 

graphical positions), /a^ yields a value 0.9932. 

Figure 7 describes the geometry of interest in calculating sunrise and 
sunset times. The half angle of darkness S at a given latitude can be 
calculated using Figure 7 as 

8 = cos ‘'“{tan 6 * tan 41 '} ( 2 ) 


d a 

where tan <J>' = — and tan 6’ = -j- . In (2), cp T represents the geocentric 
latitude of the point of interest and 6 * represents the geocentric latitude 
of the sub-solar point for a particular day (declination of the sun) . Using 
4 and 5 as geodetic latitudes, ( 1 ) can be used to rewrite ( 2 ) as 

8 = cos ^j(. 9932)^ tan 5 tan <f> } . ( 3 ) 


The longitude of the sub-solar point A g for a given day can be obtained 
from the ephemeris of the sun (see ref. 17) by using the equation of time 
tabulations . Thus 




Time in minutes 


4 min/degree of longitude 




The ephemeris also provides the declination of the sun. 
Using Fig. 7, the longitude of the sunrise line is 
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Fig. 7. Geometry of interest for calculating sunrise and sunset times. 


Sun 
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X SR = X S + * ' 3 

while the sunset line is 

x ss A s + 71 ” ® * 


Therefore, the GMT of sunrise and sunset can be calculated as 

X " X SR 

GMI SR = 15 + 12 

X ' X SS 

GKr ss = — IT- + 12 


(4) 


where the denominator on the right side of (4), 15, represents degrees per 

hour of earth rotation and the sub-solar point is defined at 1200 GMT. 

As an example to determine the GMT and GMT _ at 71°W, 42°N on June 13, 

SK So 

use (3) to get 

6 = cos‘ 1 [(.9932) 2 tan(23.16°) tan (42°)] 


3 = 67.7° . 


The subsolar time is 0.3 min after 1200 GMT, so that 

A = - ~ = -0.075° . 

S 4 

Here A CD = -0.075 + 180 - 67.7 = 112.2, 

and A ss = 247.6° 


Therefore, 


at 71°W, 


42°N, (4) 
gmi sr = 
Glfr ss = 


yields 


71 - 112.2 
15 


71 - 247.6 
15 


12 

12 


0915 
0012 . 
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4.3 Daytime and Nighttime Phase Levels 

With respect to the T-model Omega phase prediction, daytime refers to 
the period of time when the receiver and selected transmitter and the inter- 
vening path are in a daylight condition. Nighttime is that period where the 
total path is in darkness. Other times are referred to as transition periods 
which are either sunrise transitions or sunset transitions, depending upon 
whether the path is changing from nighttime to daytime or vice versa. 

The ionosphere level is higher at night than during the day so that the 
phase velocity of the Omega signal is lower at night than during the day. This 
diurnal perturbation of phase velocity changes the wavelength of the Omega signal 
so that at a given receiver position the measured nighttime phase will be greater 
than the measured daytime phase. The T-model assumes that the nighttime phase 
is constant and greater than the assumed constant phase during the daytime 
period. A linear change in phase is assumed during the transition periods. 

Figure 8 illustrates a typical T-model phase prediction curve. In determining 
the daytime and nighttime levels, measured estimates of phase velocity 
relative to the speed of light have been used. In ref. 13, relative phase 
velocities are estimated as ( c / v p ) night = 1.00040 and (c/v p > day = 0.99730 
at 10.2 kHz. In Figure 8, ^ represents the free space wavelengths from 
transmitter to receiver, and <f>^ represents the chart wavelength corresponding 
to the average phase velocity v p = c/.9974 used by the Navy in published 
Omega charts. 

The nighttime phase prediction is then calculated as 

^night 

where 

6 

r c 

Here, is the free space phase 
by a distance d, with c the free 
f the frequency of transmission. 

^day 


k) 

\ p/night 


2jTd 

A 


( ¥ ), 


at the receiver separated from the transmitter 
space phase velocity (speed of light) , and 
Similarly, the daytime phase prediction is 


= Qday 
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The sunrise transition phase prediction is given by simply drawing a straight 
line from the nighttime phase value at the earliest sunrise time (0915 GMT 
at Cambridge) to the daytime phase prediction at the latest sunrise time 
(0940 GMT at Trinidad) . The sunset transition phase prediction is obtained 
by extending a straight line from the daytime value at the earliest sunset 
time (2230 GMT at Trinidad) to the nighttime phase prediction at the latest 
sunset time (0012 GMT at Cambridge) 



Fig. 8. T-model phase prediction at Cambridge, Mass., for 10.2 kHz Omega 
signal transmitted from Trinidad transmitter on June 13. 

4.4 T-Model Phase Difference Predictions 

In predicting the phase difference or LOP type phase measurement, the 
T-model is used to provide a phase prediction at the receiver from each 
transmitter. The difference between the two curves is then used as the 
phase difference prediction at the receiver site. Figure 9 illustrates a 
typical LOP phase prediction for a receiver located at Cambridge, Massachusetts. 

The T-model provides directly a phase prediction at a given receiver site. 
In order to be of practical use, a phase correction is needed which can be 
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Fig. 9. T-model phase and phase difference prediction at Cambridge, Mass., 
10.2 kHz Omega signal transmitted from Trinidad and from Hawaii on 
June 13. 



applied to an actual received signal to reduce navigation errors. To obtain 
a sky-wave correction (SWC) , the T-model predictions are compared to the chart 
phase, which is the nominal phase for a particular path based on a nominal 
value for phase velocity. This value, v^ = c/.9974, is the velocity used 
in navigation charts (see ref. 1) published by the Navy. 

Figure 10 illustrates the sky-wave corrections for stations B and C and 
the correction for (B-C) LOP at Cambridge, based on the predictions illustrated 
in Figure 9. 

4.5 Free Space Wavelength Distance 

To calculate the daytime and nighttime phase prediction at a given 
receiver, it is necessary to determine the distance between the receiver and 
the transmitter. Distances between points on the earth's surface can be 
approximated by transforming the geodetic coordinates of the points to 
geocentric coordinates and then computing the arc length as if it were a 
geodesic on a sphere. 

Figure 11 illustrates two points P and P' on the earth's surface with 
their geocentric coordinates (a, 3) and (a’, 3*) respectively. By constructing 

two vectors from the earth's center to each of these points, the cosine of 

the angle between these vectors can be related to their coordinates by their 
dot product. Define the vector 

— A A 

P = i R cos a cos 3 + j R sin a + k R cos a sin 3 
and P' = i R cos a* cos 3 1 + j R sin a' + k R cos a' sin 3* 

Then y = cos 

- cos ^{cos a cos a' cos 3 cos 3’ + sin a sin a 1 

+ cos a cos a T sin 3 sin 3') 

The distance between the two points P and P ? , d, can then be found as 

d = Ry 

where R is the mean earth radius and y is in radians. 
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Fig. 10. T-model sky-wave correction predictions at Cambridge, Mass., for 
10.2 kHz Omega signal transmitted from Trinidad and from Hawaii 
on June 13. 
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Fig. 11. Calculating distance y on the earth's surface (between 
two points) . 


4.6 Effect of Transmitter-Receiver Distance on Phase Corrections 

As has been shown, the trapezoidal model prediction of phase measurement 
can be used to generate sky-wave corrections to be applied to phase measure- 
ments made at the receiver. This analysis relates the error in estimating the 
distance between transmitter and receiver to the error in the sky-wave 
corrections yielded from the trapezoidal model. 

4*6.1 Daytime and nighttime error .— Consider a receiver transmitter 
path. Let d^ be the distance between transmitter and receiver, where 

d i = d Tl + E dl • 

with d Tl the true transmitter receiver path difference and e dl the distance 
error in estimating d^. 

The nominal phase reading at the receiver would be 


= 0 + <f> 


M. (A + . S _ 2TTf . 2irf 

\ \ v ■ % \ * “ 
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where represents the nominal phase reading based on the trite distance, 
d T ^, and <J> ei represents the phase error introduced by the error in the distance 
estimate* The quantity v^ represents the nominal phase velocity &t the 

velocity, where v c = 0.9974c, with c representing the velocity of light. 

Using the trapezoidal model, 


+i. 

day 



where <j> c is the phase at the receiver based on a free-space wavelength. The 
nominal phase reading is 


fe) 


♦i - Mr) ? - 9974 * c 


( V p)c 


Here I — } is a daytime average relative phase velocity, 
/day 


Therefore, 


day 


.9974 Iv ) A 

\ p/day 


and similarly, 

* = /c_\ 

1 night ,9974 Vp/night 

where = 0.99730*, and p- J = 1.00040*. 

V P/day \ V p/night 

The phase correction or SWC then becomes for daytime 



*Nominal based on Pierce’s results. 
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The first term in the expression for A<J> is the phase difference based on 

day 

the true distance between receiver and transmitter, whereas the second term 
represents the contribution to phase error from the error in the distance 
estimate. Then 



where <{> T ^ is the nominal phase measurement based on the true distance 
For nighttime, 




To convert e . , to n. mi. 
A<f> 

receiver as 


, define distance error incurred at the 
X 

e = -£ e 

1 2 it e A4, 


where A^ is chart wavelength in n. mi. Therefore, 



where ^ v d T - -j— d . Recalling that e = ad , then 
lclcl a l i l 
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s = i _ Aj>24.?y B 

w L - 997 *J d i 


and similarly 


night 


Jg ffll ght _ J 

.9974 J E d 1 ' 


For receiver-transmitter pair 2, 


.9974 e d. 


'night 


pzaight _ I - 
.9974 X l e d 2 • 


Then the error incurred in the phase difference SWC generated from the 
trapezoidal model is 


a 1.2 

Ae = e_ - 
day 1, 2, 

day day 


,-P/ da 2 

.9974 


‘■*1 ■ v 


'night 


ft), 


p /night 
.9974 


The errors and are independent so that 


2 2 

o. 0 - a 
1>2 e 


+ or a. 2 = [a 2 + a 2 1 

% 1>2 K \\ 


and the standard deviation of the daytime and nighttime errors becomes 


ft), 


_ , \ p/day 2 2 

k 1 “ orn/ a + a 

1,/ day L J L d i d 2 . 


“night 


P /night 
.9974 


2 2 

cr + a 

d l d 2 
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These represent the maximum SWC error incurred in navigation due to errors 
in the estimates of the distances between transmitter and receivers when 
using the trapezoidal model to generate sky-wave corrections . 

4.6.2 Error as a function of time of day .— The section above has 
described the maximum position error using trapezoidal model SWC as a function 
of the error in estimating the distance between the receiver and transmitter 
pair. Next, consider this error as a function of time of day. 

Define 

AT _ = difference in time between sunrise at the 

bK 

receiver and transmitter 

ATgg = difference in time between sunset at the 
receiver and transmitter. 

These times are indicated in Fig. 12. 



From Fig. 12, slope m^ is given by 


m l = 


* (.-) - 4 > ( £ ) 

c V c v , 

n d 


AT 


SR 


and 


m 2 = 


4> (-) - 4 (-) 

C V c V , 

n d 


AT 


SS 


TEiT SR 


TEAT SS ‘ 
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For the night-to-day transition during the period AT 

SR 


A-KT) = 

4>(T) - ^ 

=3 

(-) + m. T - d). 

T C V 1 r l 




n 


TeAT SR 


A*(T) 


£) - * £) T 


C V 


4 > (”) - 

cV 

n 


n 


c v 


AT 


SR 


- <f>. 


TeAT 


SR 



TeAT sr 


where is chart phase and (— ) = (•£— ) and (—) = ) 

V A V V v V 7 

d p day n p night 


For T = 0; A4»(T) « <J> 


For T = AT sr ; A$(T) = <J> 


1 I .9974 M 

LJk.1 

l[ .9974 J 


Fig. 13 illustrates A0(T) as a function of TeAT c _. 

SR 
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For the day-to-night transition period in Fig. 14, 


A*’cr> = ken - = k c (^> d + v ' *1 



Eig. 14. Phase error in day-to-night transition period. 

Figure 15 presents the position error as a function of time-of-day superimposed 
on the trapezoidal prediction of Fig. 12. 
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Fig. 15. Position error as a function of time of day for receiver- 
transmitter path 1. 


Consider a transmitter pair where there is uncertainty in each of the 
receiver-transmitter distance measurements. Figure 16 represents the phase 
difference prediction using the trapezoidal method and the associated position 
error standard deviation using the transmitter pair. 

As can be seen in Fig. 16, the maximum error in position estimate using 
trapezoidal SWC as a function of receiver-transmitter distance measurements is 
during the total nighttime period. The rms error is given by 



(ID 


where 



= variance of distance error in path 1, 


= variance of distance error in path 2, and 


- standard deviation of position estimation error based on 
phase difference measurements. 
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4.6.3 Analysis .— As can be seen from this analysis, the error in position 
estimation from phase difference measurement employing trapezoidal SWC at a 
receiver is relatively insensitive to errors in distance measurement between 
the receiver and the transmitters used. (Trapezoidal SWC are thus relatively 
insensitive to error in receiver-transmitter distance measurement.) The 
maximum error occurs during nighttime conditions, where nominally the ratio of 
(c/v) n = 1.0004 so that the standard deviation of position error 



.003 a 2 




( 12 ) 


Thus, an rms error of 18 n, mi. can be sustained in estimates of receiver- 
transmitter distance and yields less than 1 cec position error at 10.2 kHz 
at the receiver. 

During daytime, an rms error of 560 n. mi. can be sustained in estimates 
of receiver-transmitter distance and yields less than 1 cec position error 
at 10.2 kHz at the receiver. 

To illustrate confidence in the method discussed in Section 4.5 for 
calculating distances between points on the earth’s surface, the distances 
from the four present Omega stations to Cambridge were computed and compared 
with values determined in ref. 13. These calculations are in terms of SWC 
chart wavelengths. 

Pierce's Geodesic 

Distance Distance 



Latitude 

Longitude 

(in Chart cec) 

Cambridge 

42°22'39"N 

(42.3775°N) 

71°07 T 03”W 

(71.1175°W) 



Norway 

66°25*15 ,, N 

(66.420833°N) 

13°09'10"E 
(346. 847222°W) 

18792.16 

18819.12 

Trinidad 

10°42 f 6.22"N 
(10. 701728°N) 

61°38 ? 20.30"W 

(61.638972°W) 

12316.34 

12311.34 

Hawaii 

21°24 , 20.67 ,, N 

(21.405742°N) 

157°49’47.75 ,, W 
(157. 829931°W) 

27747.83 

27753.84 

New York 

43°26 f 40. 92"N 
(43.4447°N) 

75°05’09.80"W 
(75 .086056°W) 

1170.86 

1171.47 


These calculations were made with the following constants: 
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one chart wavelength = P . r . ee space wavelength at 10.2 kHz 

.9974 

= 2.946641 x 10 4 meters 

one mean earth radius = 6371 km 

= 216.212271 chart wavelengths. 

The largest difference between the geodesic and the Pierce values is in 
the Norway-Cambridge distance, which is off by 26.96 cec, or 7.944 km. 
The largest (daytime) error that this can cause when computing sky-wave 
corrections using the T-model is, from (5), 


e,. . = .0001 x 26.96 « .0027 cec 

distance 

4.7 Summary 

To obtain the T-model phase correction at a particular receiver site, the 
latitudes and longitudes of the receiver and transmitter are needed. The 
ephemeris of the sun, along with an estimate of the receiver-transmitter 
distance, can then provide all the information needed to get the sky-wave 
correction at a particular frequency. Constants needed are relative daytime, 
nighttime, and nominal phase velocities and estimates of the earth radius. 

The T-model method offers a very simple scheme for generating sky-wave 
corrections. The resulting correction is constant during the receiver- 
transmitter path daytime and nighttime periods. During the transition 
periods, the correction is essentially the daytime correction modified 
according to the percentage of the path in darkness. The change in phase 
correction as a function of the percentage of the path in darkness is assumed 
linear. 

Appendix D provides a flow chart description of a computer program to 
generate phase corrections using the T-model scheme. 
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5.0 COMPARISON OF TRAPEZOIDAL MODEL, NAVY MODEL, 

AND ACTUAL DATA 

This chapter provides a comparison of predicted phase difference measure- 
ments at 10.2 kHz over a one-year period with actual phase measurements made 
by J. A. Pierce at Harvard University in Cambridge, Massachusetts. Phase 
difference predictions are provided from two sources. The published Navy 
sky-wave corrections have been used to obtain one set of predicted phase 
difference measurements by applying the published corrections to the nominal* 
phase prediction in the Cambridge area using the Trinidad-Hawaii (B-C) 
transmissions. A second set of phase difference predictions has been obtained 
using the trapezoidal model described in Chapter 3.0. Since the Navy 
predictions are provided as two-week averages, the T-model predictions used 
are 15-day average phase difference predictions. The daily phase measurement 
Pierce data have also been averaged over 15-day periods. One 15-day period 
for each month from October 1970 through September 1971 is considered. All 
periods are the first 15 days of the month except for May 1971, when the 
Pierce data were incomplete. For May 1971 the average of the last 15 days 
of Pierce data is used. 

Figures 17-25 compare phase difference predictions from the published 
Navy tables (N) , the T-model (T) , and the Pierce data (P) for the months 
Jan. -Sept. 1971. Figures 26-28 provide the same comparison for the period 
Oct. -Dec. 1970. Predictions and measurements are provided in centicycles at 
10.2 kHz averaged over 15 days for each hour of the 24-hour period. Lines 
have been drawn between points to assist in reading the plots. 

Figures 29-40 provide pairwise differences between the curves given in 
Figs. 17-28, In each of Figs. 29-40, the curve defined by points labeled "X" 
represents the difference between the Pierce average measurements and the 
T-model average predictions for each hour of the 24-hour day. Curve n Y" 
represents the difference between the Pierce average measurements and the 
Navy predictions. Curve "Z" represents the difference between the two 
predictions. Relative phase difference is indicated in centicycles of 10.2 kHz. 

Plots of the monthly mean difference between each pairwise comparison 
are given in Figs. 41-43. The monthly mean is determined by computing the 
24-hour average of the hourly values given in Figs. 29-40. 

*Nominal phase is based on chart velocity. 
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Figure 41 illustrates the mean difference between the Pierce measurements 
and the T-model predictions for each 15-day period. Figure 42 shows the mean 
differences between the Pierce measurements and the Navy predictions, while 
Fig. 43 provides the mean differences between the T-model and Navy phase 
predictions. Additionally, these three figures show three-month average 
differences for each pair to indicate seasonal variations. 

Figures 44-46 provide monthly and quarterly rms difference for each 
pairwise comparison. These rms values for each month are determined with 
respect to the monthly means given in Figs. 41-43 for pairs X, Y, and Z, 
respectively. The quarterly rms difference is provided to indicate seasonal 
variations. 


5.1 Analysis 

In Figures 17-28 it can be seen that the T-model phase predictions do 
provide a generally accurate representation of what the phase difference 
readings will be. The diurnal shifts are definitely present; however, the 
largest error seems to be in prediction of the time of day when the diurnal 
effects occur. The all-daytime and all-nighttime periods are illustrated by 
the horizontal straight lines on the trapezoidal predictions. Variations 
within these periods as evidenced by the actual data and as accounted for 
by the Navy predictions are, of course, not obtainable with the T-model 
predictions. After some examination, it appears that some simple filtering 
of the T-model predictions to smooth out the "sharp corners" might improve 
these predictions. 

Analyzing Figures 29-40, it is clearly evident that the largest errors 
in the T-model predictions occur during the transition regions. The largest 
errors are on the order of 25 cec, occurring in Jan., Feb., and Mar. during 
the nighttime-daytime transition on the Cambridge-Hawaii path. A shift in the 
calculated sunrise times at Hawaii and Cambridge could improve this error 
considerably. 

Generally , from Figure 41, the T-model predictions are too low in the 
winter months (Jan. -Feb.) and are high for the remainder of the year. 

Closest agreement occurs in Jun. and Jul. and shows very small average error 
(less than .3 cec) with the summertime predictions less than 1 cec on the 
average. Considering the mean values above, the T-model predictions do as 
well or better than the Navy predictions (Figure 42) with the yearly average 
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deviation from the measured value less than 3 cec for both. Comparison of 
the Navy and T-model predictions in Fig. 43 shows that the trapezoidal 
predictions are generally higher than the Navy predictions. 

The rms error in the T-model predictions in Fig. 44 illustrates why 
the T-model predictions have been objected to previously (see ref. 14). 

The rms error as compared to that of the Navy predictions (Fig. 45) is some 
3 cec higher comparing the yearly average rms error. This error is the error 
when comparing the predictions with the measured phase values. As illustrated 
In Fig. 46, the rms error between the two predictions averages out to 
approximately 8 cec. However, comparison of each method of phase prediction 
to the measured phase is a more valid comparison. 

5.2 Conclusions 

The T-model phase predictions account for the most significant Omega 
phase perturbations caused by diurnal ionospheric variations. The greatest 
advantage of this method is the inherent simplicity of generating the pre- 
dictions and/or corrections. As has been illustrated, the incurred first-order 
error is comparable to that achieved with the complex Navy sky-wave correction 
method. Although the second-order error of the T-model is not as good as 
the Navy predictions, it appears that there may be a possibility of filtering 
the trapezoidal predictions to improve the second-order errors. Furthermore, 
the simplicity of generating the T-model predictions enhances the possibility 
of real-time employment in the SWC gradient corrected mode of differential 
Omega (see Section 6.3). 
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Fig. 17. 10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC Tables (N) , Trapezoidal 

Model (T), and Pierce Measurements (P) for period 1-15 January 1971. 
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Fig. 18. 10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC Tables (N) , 

Trapezoidal Model (T), and Pierce Measurements (P) for period 1-15 February 1971. 
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Fig. 19. 10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) 

Trapezoidal Model (T), and Pierce measurements (P) for period 1-15 March 1971. 
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Fig. 21. 10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) , 

trapezoidal model (T), and Pierce measurements (P) for period 16-30 May 1971. 
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Fig. 22. 


10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) , 
trapezoidal model (T), and Pierce measurements (P) for period 1-15 June 1971. 
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Fig. 23. 10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N), 

trapezoidal model (T), and Pierce measurements (P) for period 1-15 July 1971. 



ThIMUAD-HAWA j 3 
M u W T *4 s B 


l.'-C 

DAYS=1-15 


15 


Ui 

LO 


-5360 

1.01 

2.0 

3.0 
<1 . 0 

5.0 
6.U 

GREENWICH 7 ■ 0 
MEAN TIME 

B.O 

9.0 
10.0 
11.0 
12.0 

13.0 

14.0 

15.0 

16.0 

17.0 
16. U 

19.0 

20.0 
21.0 
22. 0 

23.0 

24.0 


'5360 


I 


K = 


PHASE D I • F"EPf Net IN CtC-S 
’5400 -5420 -5440 -5460 -54bO 

N*. k ! I I 



-55?0 

I 


-5540 

I 


-5560 

I 


LEGEND 

= PIERCE DATA 
= TRAPEZOID CORRECTIONS 
“ NAVY SWCb 

= PIERCE AND TRAP. COINCIDENT 

- TRAP. AND SWC 

- PIERCE AND SWC 

- ALLCOINCIDENT 


Fig. 24. 10,2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) 

trapezoidal model (T) , and Pierce measurements (P) for period 1-15 August 1971. 
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Fig. 25. 


10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) , 
trapezoidal model (T), and Pierce measurements (P) for period 1-15 September 1971. 
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10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) 
trapezoidal model (T), and Pierce measurements (P) for period 1-15 October 1970. 
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Fig. 27. 10.2 kHz Omega phase predictions for B-C at . Cambridge , Mass., using Navy SWC tables (N) 

trapezoidal model (T), and Pierce measurements (P) for period 1-15 November 1970. 
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Fig. 28. 


10.2 kHz Omega phase predictions for B-C at Cambridge, Mass., using Navy SWC tables (N) , 
trapezoidal model (T), and Pierce measurements (P) for period 1-15 December 1970. 
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Fig. 29. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 January 1971. 
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Fig. 30. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 February 1971. 
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Fig. 31. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) » and trapezoidal prediction minus Navy 
prediction (2) at Cambridge, Mass., for period 1-15 March 1971. 
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Fig. 32. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (2) at Cambridge, Mass., for period 1-15 April 1971. 
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Fig, 33. B-C 10,2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X) , Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 16-30 May 1971. 
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Fig , 34. 


B-C 10 2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 

^ * ?± * rce m * asur&ment minus Nay y Prediction (Y) , and trapezoidal prediction minus Navv 
prediction (Z) at Cambridge, Mass., for period 1-15 June 1971. 
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Fig. 35. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y), and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 July 1971. 
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Fig. 36. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X) , Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 August 1971. 
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Fig. 37. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y), and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 September 1971. 
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Fig. 38. 


prediction at Cambridge, Mass., for period 1-15 October 1970. * 
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Fig. 39. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 November 1970. 
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Fig. 40. B-C 10.2 kHz Omega phase difference comparison of Pierce measurement minus trapezoidal prediction 
(X), Pierce measurement minus Navy prediction (Y) , and trapezoidal prediction minus Navy 
prediction (Z) at Cambridge, Mass., for period 1-15 December 1970. 
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Fig. 41. Monthly averages (using one 15-day period per month) B-C phase 

difference for Pierce measurements minus trapezoidal predictions 
of 10.2 kHz Omega phase at Cambridge, Mass., for Jan. -Sept. 1971 
and Oct. -Dec. 1970. 
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Fig. 42. Monthly averages (using one 15-day period per month) B-C phase 

difference for Pierce measurements minus Navy predictions of 10*2 
kHz Omega phase at Cambridge, Mass., for Jan-Sept. 1971 and Oct.- 
Dec. 1970. 
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Fig. 43. Monthly averages (using one 15-day period per month) B-C phase 

difference for trapezoidal predictions minus Navy predictions of 
10.2 kHz Omega phase at Cambridge, Mass., for Jan. -Sept, 1971 
and Oct. -Dec. 1970. 
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Fig. 44. Monthly rms (using one 15-day period per month) B-C phase difference 
for Pierce measurements minus trapezoidal predictions of 10.2 kHz 
Omega phase at Cambridge, Mass., for Jan. -Sept. 1971 and Oct. -Dec. 
1970. 
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Fig. 45. Monthly rms (using one. 15-day period per month) B-C phase difference 
for Pierce measurements minus Navy predictions of 10.2 kHz Omega 
phase at Cambridge, Mass., for Jan. -Sept. 1971 and Oct. -Dec, 1970. 
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Fig. 46. Monthly rms (using one 15-day period per month) B-C phase difference 
for trapezoidal predictions minus Navy predictions of 10.2 kHz Omega 
phase at Cambridge, Mass., for Jan. -Sept. 1971 and Oct. -Dec. 1970. 
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6.0 DIFFERENTIAL OMEGA 


6.1 Differential Omega Concept 

The concept of Differential Omega was first formally proposed in 1966 
by the Omega Implementation Committee set up by the Department of the Navy 
(see ref. 18). The assumption made is that the phase difference error of 
the Omega signals within any small region (approximately 300 n. mi. in radius) 
will be very nearly constant. The assumption that sky-wave corrections are 
constant over a 240 n. mi. square region is made in the published tables, 
where one correction applies over a 4° x 4° area. 

One mode of operation would be to locate a base receiver at a fixed 
known position so that true Omega phase is known. This receiver can then 
at any time make an actual phase difference measurement and obtain a phase 
difference measurement correction applicable to other receivers in the 
differential region. Some auxiliary means of communications could then be 
used to transmit the correction information from the base receiver to other 
receivers within the differential region. 

The method offers a real-time measurement of the phase perturbations 
within a relatively small area of operation about the base receiver. This 
should serve to reduce navigation error over the ordinary Omega navigation 
procedure where phase perturbation predictions are used. Differential Omega 
studies, including refs. 19-28, have shown that the differential method can 
achieve better navigation accuracies than ordinary Omega. 

6.2 Navy SWC Table Analysis 

To gain insight into the nature of time and space correlation of Omega 
phase differences, a simple analysis of the Navy sky-wave correction (SWC) 
tables has been carried out. Using the tables, lines of constant phase 
correction for given transmitters and lines of constant phase difference 
correction for various pairs of transmitters have been plotted in the North 
American region. Figures 47 and 48 provide plots of constant Navy SWC for 
transmitters B and D, respectively, at 0700 GMT for the 1-15 July period. 

The value of phase correction given in the SWC tables is assumed to hold at 
the center of the 4° x 4° lattice square. Linear interpolation is used between 
lattice regions in determining where these isolines are plotted. Figure 49 
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values for each 4°x4° lattice grid. 












































shows isolines of phase difference correction for the station pair B-D obtained 
from Figs. 47 and 48. Figures 51 and 52 show isolines of phase difference 
correction for station pair A-C at 0600 GMT and 1700 GME, respectively, for 
the period 1-15 July. Isolines are shown for the area immediately around 
Langley Research Center (LRC) only. From these plots it can be seen that 
changes in phase difference corrections vary with direction and are generally 
non-linear in any given direction. Further analysis has been carried out in 
the immediate area of LRC by calculating the gradient of phase difference 
correction for each hour of the day using the published SWC tables for 1-15 July. 
The gradient as used here is defined as the vector which represents the magnitude 
and direction of greatest phase difference correction change for the SWC table 
region in which LRC is located. The following illustrates the technique used. 

The set of 9 regions in Fig. 50 represents the sky-wave correction region 
containing LRC (a 22 ) and the 8 surrounding it. For each d x d region and each 
Omega transmitter, there is a published SWC. The a^’s represent the differ- 
ence between any two corrections. It is desirable to find the gradient of a 
field represented by the a^'s and centered at a 22 ; ^* e *» t ^ ie direction and 
magnitude of the largest change in the r s around a 22 » 
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Fig. 50. Localized Omega SWC grid for determining gradient of 
sky-wave corrections at center grid. 
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If we consider a surface containing points raised from the center of each 
region by an amount equal to a_ for that region, then the closest approxima- 
tion to the partial derivatives in the gradient will be for the East 
coordinate 
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Fig. 53. Magnitude (a) and direction (b) of gradient of 10.2 kHz Omega SWC 
difference for A-B phase difference measurements for area centered 
at Hampton, Va., for each hour during the period 1-15 July. 
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and for the North coordinate 


Grad 
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a 12 a 32 

2d 


These orthogonal components then determine the magnitude and direction of 
the resultant vector which is the gradient. The above approximation is 
considered to be a good one, since the gradient varies slowly with distances 
on the order of d. 

Figures 53 through 58 are plots of the magnitude and the direction of 
the gradient centered at Hampton, Virginia, for the 10.2 kHz SWC differences 
for each hour of the 24-hour day. All possible pairs for stations A, B, C, 
and D are shown. The magnitude of the gradient is shown in centicycles (cec) 
per degree of latitude or longitude, while gradient direction is given in 
degrees measured ccw from East. It can be seen that for most station pairs 
a diurnal variation is evident in the gradient magnitude and direction. 

This is most pronounced with pair B-C in Fig. 57, where the baseline between 
the Trinidad and New York transmitters is parallel to the sunrise-sunset line 
and Langley Research Center is located very near the baseline. 

This same information is presented in Figs. 59 through 64 in polar form, 
superimposed on an Omega map of part of the Eastern United States. In these 
figures, individual points indicate the termini of each hourly gradient for 
the two-week period analyzed. As mentioned previously, each hourly gradient 
vector has a magnitude and direction as indicated in Figs. 53 through 58. The 
vectors actually drawn in the figures represent median gradients for various 
periods in the day where the gradient remains relatively stationary. 
Additionally, the daily average gradient, Z, is shown. This is calculated 
by determining the average East coordinate and the average North coordinate 
for all hourly gradients over the 24-hour period. A trend is evident 
in that the average gradient is generally almost perpendicular to the station 
pair LOP at LRC for each station pair; this is not necessarily the situation 
for given individual hourly gradient directions. (Note that the Omega map 
only shows LOP’s for A-C, A-D, B-C and B-D, which are considered the most 
usable pairs in the area.) 

The magnitude of the gradients is generally less than about 2 cec per 
degree of latitude or longitude, which means that variations of phase difference 
measurement of less than 10 cec (.8 n. mi.) would be generally true within a 
300 n. mi. radius of LRC. This would allow for differential Omega navigation 
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Fig. 55. Magnitude (a) and direction (b) of gradient of 10.2 kHz Omega SWC 
difference for A-D phase difference measurements for area centered 
at Hampton, Va„, for each hour during the period 1-15 July. 
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Fig. 56. 


Magnitude (a) and direction (b) of gradient of 10.2 kHz Omega SWC 
difference for B-C phase difference measurements for area centered 
at Hampton, Va., for each hour during the period 1-15 July. 
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Fig. 58. Magnitude (a) and direction (b) of gradient of 10.2 kHz Omega SWC 
difference for C-D phase difference measurements for area centered 
at Hampton, Va., for each hour during the period 1-15 July. 
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Fig. 59. Median gradients of 10.2 kHz Omega SWC differences for A-B phase dif- 
ference measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for the 
period is shown as £. Vector plots are superimposed on Omega LOP chart. 
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Fig. 60. Median gradients of 10.2 kHz Omega SWC differences for A-C phase dif- 
ference measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for 
period is shown as £. Plots superimposed on Omega LOP chart. 
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Fig. 61. Median gradients of 10.2 kHz Omega SWC differences for A-D phase dif- 
ference measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for the 
period is shown as E. Plots superimposed on Omega LOP chart. 




Fig. 62. Median gradients of 10.2 kHz Omega SWC differences for B-C phase dif- 
ference measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for the 
period is shown as E. Plots superimposed on Omega LOP chart. 
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Fig. 63. Median gradients of 10.2 kHz Omega SWC differences for B-D phase dif- 
ference measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for the 
period is shown as Z. Plots superimposed on Omega LOP chart. 
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Fig. 64. Median gradients of 10.2 kHz Omega SWC differences for C-D Phase differ- 
ence measurements for indicated hourly periods for area centered at 
Hampton, Virginia, during the period 1-15 July. Mean gradient for the 
period is shown as E. Plots superimposed on Omega LOP chart. 
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errors of less than 1 n. mi. within this "differential region." 

Although this analysis may not provide a reliable quantitative measure 
of the navigational accuracy obtainable with differential Omega, it does 
provide some qualitative measure of the spatial and temporal correlation 
of Omega corrections within a relatively small area. As an example, for 
station pair B-D, SWC differences may be adequately accounted for by measuring 
a daytime value and a nighttime value, whereas for another pair, B-C, more 
frequent phase difference corrections would be needed to maintain good 
navigational accuracy. 


6.3 Modes of Differential Omega 

There are several modes of operation with differential Omega that have 
been conceptualized. Common to these modes is the need for a base station 
which is in a precisely-known position and which is able to determine the 
difference between the actual received phase difference and the charted 
phase difference for a given transmitter pair. 

For aircraft usage, ref. 22 has defined a differential mode and a 
sky-wave corrected differential mode. In the differential mode, the difference 
between the measured phase difference and the charted phase difference can 
be communicated from the base station receiver to all subscribers within the 
differential region. These subscribers could then directly apply the phase 
corrections to measured phase differences. In the sky-wave corrected mode, 
a larger area of operation might be possible. Here, the base station applies 
the published SWC to the measured phase difference for any given transmitter 
pair and then determines the difference between the corrected phase difference 
and the charted phase difference. This difference is then communicated to 
subscribers who would first apply the appropriate published SWC to their 
measured phase difference and then apply the base station-determined correc- 
tion. This might have application over fairly large regions to improve 
navigational accuracies, particularly during times when anomalies such as 
sudden ionospheric disturbances (SID's) occur. 

Another mode of operation involves use of differential Omega to track 
drones. (See ref. 28.) Here, the vehicle requiring position location is 
equipped with an Omega receiver and a transmitter to telemeter the phase 
measurements to a base tracking station. The tracking station then applies 
the measured Omega phase difference correction to the telemetered phase 
difference measurements and can use differential Omega to track the drone. 
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This particular mode might also have application in a Search and Rescue 
system. 

A variation of the SWC differential mode, the "sky-wave correction 
gradient corrected differential Omega," is conceivable particularly for 
large-area usage of differential Omega. In this mode, the gradient of sky- 
wave corrections at the base station receiver is used to adjust the differential 
correction at the mobile receiver according to the distance between the two 
receivers and the azimuth from base to mobile receiver. This mode is 
equivalent to the SWC differential mode; however, the actual sky-wave 
correction at either site is not required — only the difference in the 
sky-wave correction at the two receivers is needed. 

Various telemetry schemes for communication between base and mobile 
receivers have been discussed and experimented with, including VHF and UHF 
radio systems with possible use of satellite relays. Appendix A includes 
some reports pertaining particularly to studies involving application of the 
differential Omega technique. 
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APPENDIX B 

THE RELATIONSHIP BETWEEN CENTICYCLES AND PHASE 

Given position T of Transmitter and position R of Receiver, then 

D = |T - Rj is the distance between the transmitter and receiver. Assume 

v p = phase velocity of signal. Then time for point of constant phase to "travel" 

from transmitter to receiver is t m = D/v secs. 

T P 

For a given frequency f a period 1/f = T can be used to define a 

c c c 

centicycle, («01) T C * Then the time t T can be expressed in terms of centicycles 
(cec) as 

t _ D/V P _ 100 D 
T cec “ <* 01 > T c ~ v p T c 


cec 


100 D 

v T 
P c 


The "travel" time relates to the phase of the received signal. A cec 
can be thought of as 2ir/100 radians so that for a given path length a given 
phase measurement can be directly obtained from the time measurement. For a 
given transmission path it was shown above that the "travel" time is defined in 
terms of the path length, the phase velocity, and the period of the transmitted 
signal. 

Assume that t^ represents the time measure of a point of zero phase 
at the receiver recorded in cec. The difference in time t M - t^ can be 
related to a phase perturbation <f> by considering two situations. If t M >t T 
then At = (t M -t T > mod 100, and 

If At<50 then <t> = - At (~) , - 180° <■ <f> < 0 and a phase lag 

condition is said to exist. 

If At>50 then <p = (100-At) (—) , 0<A <p <_ 180° and a phase lead 

condition is assumed. 

For W At = <VV mod 100 > and 

If At>50 then 0 = (At-100) (“•) , - 180° < (| < 0 and a phase 

lag condition exists. 

If At <50 then <f> = At (^-) , 0 < <t> < 180° and a phase lead 

condition is assumed. 

For this situation the term "lane width" is the maximum distance over 
which the phase angle can be determined unambiguously. For this time scheme 
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variations of phase on the range [-180° <_ $ £ 180°] can be determined. 
Obviously, a phase lead of 200° cannot be differentiated from a phase 
lag of 160° so that 360° or 1 cycle is the actual lane width. 

Table B-l shows the lane widths W^(f) as calculated, assuming that 
the phase velocity corresponds to the speed of light c, for frequencies 
of interest in Omega. The equation 


V f > 


V f > 

f 


provides the lane width as a function of phase velocity and frequency, f. 


Table B-l. Omega frequencies, frequency differences and 
associated periods. 


Frequency 

Velocity of 
Constant Phase 

Lane Width 
Nautical Miles 

Period T (f ) 

1.1| kHz 

161783 n. mi. /sec 
c 

142 n.mi. 

.88 msec 

3.4 kHz 

c 

48 n.mi. 

.294 msec 

10.2 kHz ' 

c 

15.8 n.mi. 

98 ysec 

11.3± kHz 

c 

14 n.mi. 

88 psec 

13.6 kHz 

c 

11.9 n.mi. 

73.5 ysec 


It should be noted here that the term "lane width" is more commonly 
used in terms of phase difference measurements between two stations along 
the baseline. Commonly "lane width" is the maximum unambiguous range 
obtained from phase difference measurements and corresponds to one-half wave- 
length along the baseline between stations. These lane widths W 1 (f ) would 

36 

be given as 
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V 


V( f > 


= = 
2f 


Vf) 

2 


and are just one-half the distances given in Table B-l. 

Off the baseline, the lane width varies in accordance with the 
hyperbolic geometry. 


121 



APPENDIX C 

TRANSMISSION TIME FOR 3.4 KHz BEAT SIGNAL 


For a waveguide propagation is characterized by a propagation constant y 
where E^, is of the form 


E 


T 


E e 


-yz 


where z is the along axis direction or direction of propagation, y is generally 
complex and is expressed as y = ct + j3, where ct is an attenuation constant and 
3* is a phase constant at a particular frequency. Phase velocity in a waveguide 
is the velocity at which a point of constant phase moves so that v^ = ^ 
represents the phase velocity. Any information signal can be thought of as 
traveling at the group velocity of the guide 


Am 



which describes the phase distortion introduced by the guide at a particular 
frequency (ref. 29). 

Consider the situation with the Omega signal where frequencies 10.2 kHz and 
13.6 kHz are transmitted. Each has a phase velocity so that 


Thus 


6 


10.2 


W 10.2 

and 


“13.6 

6 10.2 

V 

3 13.6 


p 13.6 

“l0.2 

and 

8 , „ . 

“13.6 


P13.6 


13.6 


P13.6 


If one were to consider modulating the 10.2 kHz signal with a 3.4 kHz signal, 
the result would be similar to the situation where one considers propagation of 
the difference frequency Af = g “ f^Q 2 t * ie actual Omega signal. A 
modulation signal, or in this case the difference frequency, can be thought of 
as information which is propagated in the waveguide. The group velocity of a 

*8 is the phase constant for a guide and usually is expressed in terms of 
radians per unit length at a given frequency. 
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• waveguide can be defined as 


Aoj 

A3 


(see refs. 30-32) 


Then 


27T(f 13.6 “ f 10.2 ) 


*3.4 


13.6 




10.2 


13.6 


10.2 


Note that D = VT so that 


D 


27T(f 


13.6 


I 


f 10.2^ 


‘3.4 2 7 f 13.6 T 13.6 


f T 

r 10. 2 i 10. 2 


) 


or 


= f 13.6 T 13.6 £ 10.2 T 10.2 

3 * 4 f 13.6 ‘ f 10.2 


Since f^ ^ ~ f -^0 2 = 3.4 kHz, and 


f 13.6 

H » AC J » 


Af 


Af 


then 


X 3.4 “ 4T 13.6 " 3 T 10.2 


which is the expression obtained in refs. 11 and 13. 

Consider reception of and f^ 2 * Using the same reasoning as previously, 

3 


T l.l| = 10 T ll| “ 9 T 10.2 • 
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APPENDIX D 

PROGRAM TO COMPARE TRAPEZOIDAL, NAVY AND PIERCE DATA 
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CALCULATION OF SUNRISE AND SUNSET TIMES 


ONE YEAR SUN 
EPHEMERIS 
EVERY FOURTH 
DAY 


LATITUDE AND 
LONGITUDE OF 
RECEIVER 



T. min 

V A sj ° - rgg3e5 ■ J ■ '-** 



LATITUOE AND 
LONGITUDE OF 
OMEGA TRANS- 
MITTERS 

*1* A 1 1 - ’> 

1 

4 


CALCULATE HALF-ANGLE OF 
RANGE OF LONGITUDE IN 
DARKNESS (WORK IN 
GEOCENTRIC COORDINATES) 


NOTE: ALL CALCULATIONS INVOLVING TINE HERE 
DONE IN MODULO 24 BUT THE LOGIC IS NOT 
PRESENTED HERE FOR SIM>LICITY. 


Bfj ” Cds' ] {(. 9932) 2 tan 6^ tan 
1-1.5 
j - 1,96 


FUNCTION SUBROUTINES 


SUNRISE AND 
SUNSET TINES 
IN GMT 


x. -x r , - v + 6. 


J 1 - 1.5 

®1T5S. j - — ^ ♦ 12 J_1 * 96 


DIFF (A.B) FINDS 


AVE (A.B.C.D) 

CALCULATE 

LEAST DIFFERENCE 

• 4 $ — 

FINDS AVERAGE 


OF 2 HOOULO 24 


OF 4 MODULO 24 

FOR EACH 

VARIABLES 


VARIABLES 



4k+l 

£ GMTSR 1 


J - (k-1 )4+l 
4k+l 

r ik = i Z 6MTSS ij 

j - (k-l )4+T 
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CALCULATION OF TRAPEZOIDAL PHASE PREDICTIONS 
FOR EACH TRANSMITTER TO RECEIVER PATH 


STORED DATA 

LATITUDE AND 
LONGITUDE OF 
TRANSMITTERS AND 
THE RECEIVER 



w 1 

*s’h = 

CALCU 
DISTAJ 
EACH 
TO TH 

■ATE 

iCES FROM 
rRANSMITTER 
RECEIVER 


TIME FOR EACH HALF MONTH, 
FOR EACH TRANSMITTER, AND 
FOR THE RECEIVER LOCATION 



Calculate 
NIGHTTIME AND 
DAYTIME PHASE 
FOR EACH PATH 
IN CEC 


where 


v rh =■ chart 
c " phase 
velocity, 

Jav = d W time 
^ phase 

vel oci ty , 

nighttime 
phase 
velocity 


night 


SET TIME EQUAL 
TO ZERO 


TYPE 1 


TYPE 4 


TYPE 2 


TYPE 3 

Tq ■ GMTSSfj 


T fl - GMTSS 1k 


T q - GMTSS 5k 


T n » GMTSS fi( , 

T 1 " 


T, - GMTSS 5k 


T, - GMTSS 1k 


Tj = GMTSS ik 

t 2 - GUTS** 


T 2 - GHTSR 1k 


T 2 = 


T 2 " GHTSR ik 

T 3 - GMTSR 1k 


T 3 = «*■* 

4 


T 3 = GMTSR ik 


t 3 3 ®n»Sk 


I 




I 


i 


ADJUST TIMES TO 
STANDARD FORM 


V T t * T 0’ 1 = M 


CALCULATE 
TRAPEZOID 
SLOPES 
SI AND S2 


T 

® 


SI 4 = 


PHIN < - PHID 1 


S2 4 = 


'1 

PHINj - PHID, 
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t = TIME 
PHI t1 = PHASE 

t - half hour 
phase values 
1,4B 

f - M 
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CALCULATE PHASE DIFFERENCES OF ALL POSSIBLE 
COMBINATIONS OF TWO TRANSMITTERS 
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APPENDIX E 


PROGRAM TO CALCULATE TRAPEZOIDAL MODEL 
FOR A GIVEN POINT ON THE 


SKY-WAVE CORRECTIONS 
EARTH 
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